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ABSTRACT 


This  thesis  explores  the  potential  advantage  of  deep  ocean  trench  placement  of  an 
acoustic  sensor  network.  The  hypothesis  is  that  a  hydrophone  deployed  in  a  deep  ocean 
trench  will  operate  in  an  environment  of  reduced  ambient  noise  such  that  passive 
maritime  surveillance  may  be  performed  with  relatively  high  signal-to-noise  ratio  (SNR). 
Barrier  or  “tripwire”  coverage  along  the  oceanic  trench  is  made  possible  by  virtue  of 
networked  acoustic  communication  among  widely  spaced  hydrophones.  Physics-based 
hydro-acoustic  propagation  modeling  supports  analysis  of  target-to-sensor  propagation 
and  networked  acoustic  communication  links  in  representative  trench  environments. 
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I.  INTRODUCTION 


Maritime  surveillance  involves  constant  vigilance  over  vast  areas  of 
responsibility.  Manned  platforms,  including  aircraft,  ships,  and  submarines,  are  scarce 
resources  requiring  costly  manpower.  They  alone  cannot  meet  the  requirements  for 
maritime  domain  awareness. 

Unmanned  systems  are  the  key  to  mitigating  this  shortfall.  In  situ  sensors  can 
provide  continuous  coverage  for  detecting  and  classifying  maritime  vessels  of  interest. 
Such  a  pervasive  and  persistent  sensing  capability,  coupled  with  near-real-time 
communication,  would  increase  the  effectiveness  of  our  manned  platforms.  For  cost- 
effectiveness,  the  autonomous  system  can  alert  on  targets  when  they  enter  a  specific  area 
or  choke  point.  Such  an  approach  allows  a  sensor  network  to  strategically  monitor  a 
larger  area  with  significantly  fewer  resources  for  sensing  and  alerting  on  targets  of 
interest. 

This  thesis  examines  the  potential  for  using  naturally  occurring  deep  ocean 
trenches  and  ridges  as  advantageous  sites  for  passive  acoustic  maritime  surveillance. 
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II.  MARITIME  SURVEILLANCE 


A.  THE  NEED  FOR  MARITIME  DOMAIN  AWARENESS 

The  safety  and  economy  of  the  United  States  depend  upon  the  security  of  the 
world’s  oceans,  thus  making  maritime  domain  awareness  (MDA)  a  vital  national  interest. 
The  continued  growth  of  legitimate  international  commerce  in  the  maritime  domain  has 
been  accompanied  by  a  corresponding  increase  in  criminal  activity.  Conflicts  are 
increasingly  characterized  by  a  combination  of  traditional  and  irregular  tactics, 
decentralized  planning  and  execution,  and  non-state  actors  using  both  simple  and 
sophisticated  technologies  in  innovative  ways  [1].  The  smuggling  of  people,  drugs, 
weapons,  and  other  contraband,  as  well  as  piracy  and  terrorism  against  vessels  and  ports, 
pose  threats  to  maritime  security  [2]. 

1.  Smuggling 

Just  as  the  world’s  oceans  are  avenues  for  international  commerce,  they  are  also 
the  highways  for  the  movement  of  illegal  commodities.  Maritime  drug  trafficking 
generates  vast  amounts  of  money  for  international  organized  crime  syndicates  and 
terrorist  organizations.  In  addition  to  using  legitimate  commercial  shipping,  traffickers 
use  special  vessels  for  transportation  of  contraband  such  as  high-speed  watercraft  and 
self-propelled  semi-submersibles  (SPSS). 

2.  Piracy  and  Terrorism 

Piracy  and  incidents  of  maritime  crime  tend  to  be  concentrated  in  areas  of  heavy 
commercial  maritime  activity,  especially  where  there  is  significant  political  and 
economic  instability,  or  in  regions  with  little  or  no  maritime  law-enforcement  capacity. 
Today’s  pirates  and  criminals  are  usually  well  organized  and  well  equipped  with 
advanced  communications,  weapons,  and  high-speed  craft  [2]. 
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Modern  terrorists  increase  their  effectiveness  and  reach  by  establishing  links  with 
sympathetic  organizations  around  the  globe.  These  terrorist  groups  use  shipping  as 
conveyance  for  positioning  their  agents,  logistical  support,  and  generating  revenue  [2]. 

3.  Law  of  the  Sea 

In  1982,  the  United  Nations  set  out  to  define  territorial  waters  and  boundaries. 
The  purpose  is  for  controlling  and  policing  of  waterways  to  prevent  illegal  fishing  and 
dumping.  The  treaty  came  about  due  to  an  increase  in  of  illegal  activities  and 
exploitation  of  resources  in  the  latter  half  of  the  20th  century  [3].  Although  the  United 
States  has  not  ratified  this  agreement,  other  nations  are  aggressively  staking  claim  to 
maritime  territory. 

4.  Anti-Submarine  Warfare  (ASW) 

In  addition  to  controlling  and  preventing  criminal  activities,  MDA  involves 
maintaining  military  dominance  to  keep  the  waterways  open  and  prevent  attacks  on  the 
homeland  [2].  A  current  ASW  thrust  is  to  develop  and  deploy  sensors  capable  of 
detecting  and  classifying  modem  diesel  electric  submarines  operating  in  littoral  waters  or 
deep  ocean  environments  [4], 

5.  Environmental  Sensing 

MDA  also  involves  monitoring  the  undersea  environment  for  purposes  of 
scientific  study  including  marine  mammal  tracking,  pollution  control,  oceanography,  and 
detection  of  seismic  events. 

6.  Interdiction 

Effective  interdiction  of  criminal  activities  and  rapid  response  to  other  threats  are 
crucial  to  the  MDA  strategy.  This  requires  a  coordinated  effort  to  share  information  and 
resources  among  military,  customs,  and  law  enforcement  communities  along  with  the 
private  sector  [2].  These  actions  facilitate  more  efficient  and  effective  interdiction 
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operations  that  are  increasingly  conducted  by  long-range,  extended-endurance  unmanned 
platforms  [5].  Success  hinges  on  the  availability  of  persistent  and  pervasive  sensors, 
along  with  near-real-time  communication. 

B.  PREVALENCE  OF  DEEP  OCEAN  TRENCHES 


Deep  ocean  trenches  exist  in  every  ocean  around  the  globe.  They  are  a  result  of 
geologic  shifts  in  the  earth’s  crust  or  are  caused  by  centuries  of  erosion  [6].  When  two 
tectonic  plates  collide,  the  younger  of  the  two  plates,  because  it  is  less  dense,  will  ride 
over  the  edge  of  the  older,  heavier  plate.  The  heavier  plate  plunges  steeply  through  the 
asthenosphere,  fonning  a  trench  on  the  upper  crust  that  can  be  as  much  as  75  miles  wide, 
more  than  a  thousand  miles  long,  and  several  miles  deep  [6]. 
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Figure  1.  Subduction  zone 


South  American  Plate 
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onvergence 

^a  with  related  features  (from  [7]). 


The  depth  of  the  trench  depends  on  the  region  and  characteristics  of  the 
subduction  zone  fault.  The  deepest  trench,  the  Marianas,  is  over  10,000  m  at  its  lowest 
point  while  other  trenches  have  depths  of  4000  m  -  8000  m  [8].  Not  all  oceanic  trenches 
coincide  with  subduction  zone  faults.  Large  height  variations  between  two  adjacent 

plates  can  result  in  instabilities  and  erosion  leaving  deep  crevasses  in  the  bathymetry. 
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Continental  plates  tend  to  be  less  dense  than  ocean  floors  and,  as  a  result,  trenches  exist 
along  coastal  boundaries  and  are  usually  much  deeper  than  the  adjacent  seafloor. 


Figure  2.  Puerto  Rico  Trench  with  perspective  view  of  the  Atlantic  Ocean  and 
Caribbean  Sea  floor,  which  shows  the  subduction  zone  cutting  along  the  north 
shore  of  Puerto  Rico,  Hispaniola,  and  Cuba.  The  surrounding  topography 
consists  of  various  undersea  canyons  and  islands,  including  the  continental 

ridge  (from  [9]). 
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1.  Surveillance  Tripwire 

Traditional  maritime  surveillance  requires  assigning  a  set  of  assets  to  an  area  of 
responsibility  (AOR).  Acoustic  sensors  can  monitor  the  ocean  for  purposes  of  detecting 
and  tracking  vessels  of  interest  while  gathering  intelligence  on  their  acoustic  signatures 
and  behaviors.  This  can  be  effective,  even  for  detecting  a  target  with  minimal 
intelligence  prior  to  interception  or  if  the  track  of  the  contact  is  sporadic.  The  downside 
of  this  approach  is  that  complete  vigilance  requires  populating  the  entire  area  with  an 
adequate  density  of  sensors  or  other  assets.  Targets  will  adopt  countenneasures  to  avoid 
detection  or  complicate  interdiction  by  law  enforcement.  Inevitably,  the  required  sensor 
density  must  increase  with  the  effectiveness  of  the  countermeasure. 

In  contrast,  a  tripwire  system  does  not  attempt  to  provide  total  coverage  for 
targets  of  interest  but  rather  will  detect  their  presence  at  a  specific  location  and  alert  the 
response  agency.  Cueing  in  this  manner  efficiently  transfers  responsibility  for 
interdiction  to  a  capable  platfonn. 
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A  major  advantage  of  this  approach  is  that  fewer  resources  are  needed  to  monitor 
strategically  significant  areas  with  the  ability  to  detect  and  classify  targets  of  interest. 
The  tripwire  approach  does  have  the  disadvantage  of  not  providing  continuous  coverage 
throughout  the  AOR.  Using  an  autonomous  sensor  network  as  a  tripwire  requires  follow- 
on  prosecution  by  other  platforms  and  personnel. 

2.  Trench  Exploitation  for  Improved  Detection  and  Noise  Reduction 

The  proximity  of  ocean  trenches  to  land  makes  them  suitable  sites  for  early 
warning  tripwires  to  detect  a  target  of  interest  moving  towards  or  away  from  shore. 
Coastal  waters  have  a  logistical  advantage  in  that  surveillance  can  center  on  areas  in 
which  targets  will  most  likely  operate.  As  Figure  4  shows,  the  required  area  coverage 
increases  with  the  radial  distance  from  port.  Therefore,  placing  the  system  in  close 
proximity  to  land  can  help  limit  the  number  of  sensors  required,  with  the  added  benefit  of 
ready  access  for  maintenance.  Coastal  waters  also  confine  movement,  due  to  landmasses 
that  remove  essentially  half  of  the  area  where  the  vessel  may  operate,  helping  to 
concentrate  the  search  area. 


Figure  4.  Area  of  required  coverage  for  adequate  surveillance  vs.  distance  from 
departure  point  (image  and  graphic  produced  using  Google  Earth). 
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Placing  the  sensor  network  close  to  land  also  has  disadvantages.  The  shallower 
waters  along  the  continental  shelf  have  decreased  detection  ranges  and  cannot  take 
advantage  of  the  deep  sound  channel  (DSC).  In  addition,  the  high  traffic  density 
significantly  increases  the  noise  floor  and  is  detrimental  to  the  target  SNR.  Multiple 
contacts  in  the  search  area  can  also  make  it  difficult  to  identify  a  singular  vessel.  Fishing 
and  dredging  in  the  near-shore  waters  can  also  disrupt  and  damage  equipment  for  the 
system. 

Subduction  zone  trenches  run  along  certain  continental  boundaries  placing  deep 
ocean  environments  close  to  land.  The  greater  depth  of  the  trench  can  be  exploited  by  a 
sensor  network  to  gain  a  large  field  of  view  for  detecting  targets.  It  is  believed  that 
placing  sensors  at  the  bottom  of  a  trench  will  reduce  noise  arriving  from  distant  sources 
through  blockage  by  the  trench  walls  surrounding  the  sensors.  This  allows  utilization  of 
acoustic  sensors  at  areas  that  are  close  to  land  while  at  the  same  time  minimizing 
environmental  noise  and  disruption  from  local  traffic. 

3.  Choke  Points 

Another  method  that  can  work,  in  conjunction  with  a  deep  trench  to  maximize 
detection  probability,  is  placing  the  sensor  network  at  an  area  of  restricted  movement. 
Tracking  a  target  can  be  facilitated  by  using  geographic  features,  such  as  landmasses, 
which  constrain  movement  and  minimize  the  possible  area  of  operation  by  the  target. 
Such  sites  are  commonly  referred  to  as  choke  points. 

C.  MARITIME  SENSOR  NETWORKS 

1.  Seaweb  Acoustic  Network 

Seaweb  is  an  acoustic  networking  technology  for  underwater  communication 
among  any  number  of  fixed  sensor  nodes,  repeater  nodes,  and  gateway  nodes.  It  is  also 
configurable  to  communicate  with  mobile  nodes  including  submarines  and  unmanned 
underwater  vehicles  (UUVs).  It  enables  an  autonomous  underwater  sensor  network  for 
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tracking  and  detecting  maritime  targets  in  an  area  of  interest.  The  challenge  in  the 
littorals  is  to  overcome  shallow  and  noisy  environments  that  can  limit  sensor  coverage 
and  communication  range  [11]. 

Seaweb’s  capabilities  depend  upon  several  factors  in  the  environment  in  which  it 
operates.  The  range  of  detection  is  primarily  related  to  the  depth  of  the  sensor.  Seaweb 
uses  bottom-mounted,  upward-looking  hydrophones  that  detect  targets  in  the  region  of 
direct-path  propagation. 


2.  Deep  Seaweb 


Deep  Seaweb  supports  the  greater  coverage  by  utilizing  the  natural  advantages  of 
the  DSC  and  the  favorable  propagation  in  open  ocean  areas.  A  deep  trench  that  runs 
adjacent  to  coastal  areas  provides  opportunities  for  utilizing  Deep  Seaweb  in  and  around 
a  shallow  area  to  improve  coverage  and  detection  as  a  surveillance  tripwire. 


Moored 

gateway 

buoy 


Surface  contact 


Biological 

contact 


Gatekeeper 

usv 


Sub-surface 

repeater 


node 

Deep  sound  channel  axis 


Tsunami  pressure  sensor  Acoustic  sensor 


Figure  5.  Deep  Seaweb  Acoustic  Network  System  Concept  (from  [12]). 
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Ill  PASSIVE  ACOUSTIC  SONAR  DETECTION 


A.  SOUND  IN  THE  OCEAN 

A  variety  of  sound  sources  contribute  to  the  ambient  noise  in  the  ocean 
environment.  Noise  from  shipping  traffic  can  travel  distances  of  1000  miles  or  more, 
depending  on  conditions  [13].  The  frequency  range  where  this  man-made  noise  is  most 
dominant  is  from  10  Hz  to  300  Hz  [13].  Noise  level  depends  on  the  shipping  density 
surrounding  a  given  area  with  close  proximity  to  shipping  lanes  and  harbors,  implying 
higher  noise  levels.  Turbulence  and  flow  past  the  sensor  can  be  a  significant  factor  in 
ambient  noise  levels  below  100  Hz  [13].  A  source  of  very  low  frequency  noise  is  the 
earth’s  seismic  activity,  which  is  not  significant  above  50  Hz  [13]. 


Figure  6.  A  composite  of  ambient  noise  spectra  (from  [13]). 
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B.  SONAR  EQUATION 


Each  term  of  the  passive  sonar  equation  contributes  to  the  Signal-to-Noise  Ratio 

(SNR). 


SNR  =  SL-NL-TL  +  DI  (1) 

Source  Level  ( SL )  is  the  energy  the  source  puts  in  the  water.  Noise  Level  (NL)  is 
the  additional  received  energy  from  external  activity  in  the  environment  as  described 
above.  NL  degrades  SNR  and  makes  the  desired  source  harder  to  detect  [14]. 

Transmission  Loss  ( TL )  is  loss  suffered  through  geometric  spreading  of  the  source 
energy  and  attenuation  in  the  medium.  The  TL  depends  on  water  column,  water  depth, 
and  sea  floor  characteristics  [14]. 

Directivity  Index  (DI)  is  the  ability  of  the  transducer  to  physically  or 
electronically  direct  its  received  response  in  space.  An  omni-directional  transducer  can 
detect  a  signal  in  all  directions  and  has  DI  =  0.  A  directional  transducer  detects  better  in 
a  certain  bearing  and  has  DI>  0  dB  [14]. 

For  a  given  sonar  system  in  a  given  noise  environment,  if  the  SNR  exceeds  the 
detection  threshold,  the  signal  can  be  detected.  The  use  of  an  acoustic  model,  such  as 
Bellhop,  can  identify  sensitivities  caused  by  changes  in  environmental  factors  [14]. 

C.  ATTENUATION 

Signal  attenuation  in  the  ocean  is  caused  by  numerous  effects  related  to  the 
chemical  makeup  of  seawater.  The  amount  of  attenuation  experienced  is  dependent  on 
the  frequency  of  the  signal,  as  seen  in  Figure  7.  Higher  frequency  energy  is  attenuated  to 
the  surrounding  environment  or  is  subject  to  certain  mechanisms  such  as  chemical 
relaxation  and  viscosity  losses  that  are  not  as  significant  at  lower  frequencies.  Equation  2 
is  a  derivation  of  the  curve  from  Figure  7  (from  [15]). 
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Figure  7.  Attenuation  vs.  Frequency  in  Seawater  (from  [15]). 
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D.  SOUND  SPEED  PROFILE  (SSP) 


Sound  speed  in  the  ocean  varies  as  a  function  of  temperature,  salinity,  and 
pressure.  These  factors  affect  sound  propagation  in  the  ocean  because  acoustic  waves 
refract  toward  slower  sound  speed  areas.  Temperature  in  an  area  can  vary  with  depth, 
location,  and  season,  but  tends  to  decrease  rapidly  with  depth  in  the  surface  layer  before 
becoming  relatively  constant.  Therefore,  sound  speed  is  generally  temperature- 
dependent  to  around  a  depth  of  1000  m  and  then  is  pressure  dependent  as  depth  continues 
to  increase.  The  resultant  shape  produces  a  slow-speed  propagation  duct  around  1000  m, 
known  as  the  DSC.  The  tendency  of  sound  to  want  to  move  toward  this  duct  results  in 
long  travel  distances  without  transmission  losses  from  vertical  spreading  or  surface 
reflection. 

Figure  8  shows  a  representative  deep  ocean  SSP  [16]  with  the  DSC  at  1000  m 
depth  and  the  sound  speed  increasing  linearly  at  depths  below  the  channel.  This  SSP  is 
assumed  in  the  modeling  of  acoustic  propagation  presented  in  subsequent  chapters  of  this 
thesis. 

The  SSP  in  Figure  8  does  not  display  a  mixed  layer  near  the  sea  surface  as  would 
likely  be  encountered  in  open  ocean  environments.  In  a  deep  ocean  operating 
environment,  where  the  sensor  is  bottom-mounted  looking  upward,  the  mixed  layer  will 
not  have  significant  effect  on  the  sensor  perfonnance  [14]. 
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Figure  8.  Munk  Sound  Speed  Profde  (from  [16]). 


E.  BELLHOP  UNDERWATER  ACOUSTIC  PROPOGATION  MODEL 

Sound  propagation  is  modeled  using  the  Bellhop  ray  tracing  and  beam  modeling 
program,  within  the  AcTUP  v2.2/“  toolbox,  to  simulate  acoustic  transmission  from  the 
source  to  the  receivers  on  the  trench  bottom  and  also  to  demonstrate  modem  operation 
within  the  sound  channel  [17].  Bellhop  is  designed  to  perform  two-dimensional  acoustic 
ray  tracing  and  beam  forming  for  a  given  sound  speed  field  with  variable  bathymetry  and 
absorption  boundaries  [18]. 
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IV.  DEEP  OCEAN  SENSING 


This  thesis  emphasizes  sound  characteristics  at  low  frequencies  where  attenuation 
is  less  detrimental  to  the  propagation  of  sound  from  the  target  to  the  sensor.  It  is  also  at 
these  frequencies  where  the  presence  of  machinery  tonals  and  propulsion  noise  can  serve 
to  classify  targets  of  interest.  Higher  frequency  sound  is  separately  considered  for 
purposes  of  acoustic  communications. 

A.  RECIPROCITY 

Modeling  of  acoustic  propagation  in  this  thesis  uses  the  reciprocity  rule  to  more 
clearly  visualize  the  characteristics  of  the  environment.  The  sensor  is  set  up  as  the  source 
in  the  modeling,  with  ray  data  signifying  paths  the  signal  travels  to  the  sensor. 

B.  RELIABLE  ACOUSTIC  PATH  (RAP) 

In  tropical  and  temperate  deep  ocean  locales,  passive  acoustic  sensing  is  most 
advantageously  performed  near  the  ocean  floor  [14].  Direct-path  coverage  in  deep  ocean 
environments  varies  with  depth  and  environmental  factors,  and  this  coverage  area 
corresponds  to  the  so-called  RAP. 

1.  RAP  Range  vs.  Ocean  Depth 

In  4000  m  water,  the  direct-path  coverage  extends  to  about  28  km  from  the 
source,  which  corresponds  to  a  7: 1  range-to-depth  ratio,  as  observed  in  Figure  9.  In  8000 
m  water,  the  range-to-depth  ratio  decreases  to  around  5:1  with  a  range  of  40  km,  as 
shown  in  Figure  10.  The  SSP  has  an  upward  gradient  below  the  DSC  axis.  As  depth 
below  the  DSC  increases,  so  does  the  vertical  component  of  the  sound  rays,  resulting  in  a 
greater  change  in  depth  for  the  same  horizontal  distance  traveled. 
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BELLHOP-  Munk  Profile  Gravel  Source  Near  Bottom  -8.  f=100Hz 


Figure  9.  Maximum  100  Hz  RAP  for  4000  m. 


BELLHOP-  Munk  Profile  -89  0.  f=10000Hz 


Figure  10.  Maximum  100  Hz  RAP  for  8000  m. 
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Detection  range  for  the  sensor  decreases  as  its  depth  approaches  the  DSC  until  the 
beam  becomes  trapped  within  the  channel.  At  this  point  the  range  becomes  limited  not 
by  direct  path  but  by  transmission  loss  from  spreading  and  attenuation. 

2.  Transmission  Loss  Effect  on  RAP 

Figure  1 1  shows  that  as  range  from  the  sensor  increases,  so  does  TL  even  at  low 
frequencies  where  attenuation  is  minimal.  Spherical  spreading  significantly  reduces  SNR 
with  an  8  dB  loss  in  the  first  12  km  of  propagation.  Depending  on  the  noise  level  in  the 
area  of  operation  the  detection  probability  will  be  limited  at  greater  ranges. 

In  deeper  water  where  the  direct-path  range  is  large,  spherical  spreading  again 
leads  to  significant  transmission  loss.  However,  deep  ocean  environments  tend  to  have 
the  advantage  of  lower  noise  levels  that  offset  this  reduced  signal. 


BELLHOP-  Munk  Profile  Bottom  Source  -89  0:  f=100Hz,  Zs=4390m,  Zr=10m 


Figure  11.  Transmission  Loss  for  100  Hz;  source  at  10  m;  receiver  at  4390  m. 
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BELLHOP-  Munk  Profile  Bottom  Source  -89  0.  f=100Hz,  Zs=3990m. 
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Figure  12.  Transmission  Loss  for  100  Hz;  receiver  at  3990  m. 


C.  TRENCH  ENVIRONMENTAL  FACTORS  AND  SOUND  PROPAGATION 

1.  Bathymetry 

Trenches  provide  conditions  that  can  be  both  advantageous  and  detrimental  to 
acoustic  exploitation.  The  depth  provides  for  greater  area  coverage  and  available  DSC 
exploitation  near  coastal  and  shallow  areas.  However,  since  subduction  zone  trenches  are 
the  result  of  recent  geological  and  seismic  events,  the  seafloor  is  uneven  and  at  times 
unstable.  Also,  bathymetric  features  may  exist  that  can  hinder  operations  by  blocking  or 
limiting  acoustic  communications. 

2.  Sound  Modeling  for  Representative  Trench  Environments 

The  shape,  size,  and  geologic  makeup  of  trenches  vary  significantly  depending  on 
location  and  no  singular  model  can  accurately  represent  the  every  environment. 
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Nevertheless,  modeling  can  generally  show  how  bathymetry  influences  the  sound 
reaching  the  sensor,  which  is  useful  for  understanding  how  the  environment  may  alter 
surveillance  coverage. 

For  the  testing  in  this  thesis,  two  trench  geometries  are  modeled  that  correspond 
to  basic  bathymetric  conditions  likely  to  be  encountered  in  a  Seaweb  deployment. 
Trench  A  is  a  symmetric  canyon  that  may  be  an  offshoot  of  a  subduction  ridge  or 
undersea  mountain  range.  The  parameters  are  26  km  for  edge-to-edge,  with  a  depth 
change  of  2330  m  and  max  depth  of  4400  m,  as  shown  in  Figure  13. 

Trench  B  (Figure  14)  is  more  representative  of  a  subduction  zone  with  a  high 
continental  wall  (8850  m  from  edge  to  bottom)  and  a  lower  open  ocean  wall  (3350  m). 
The  cross-sectional  width  is  1 13  km  from  edge  to  edge. 

Both  trenches  are  characterized  by  reflective  rock  walls  that  have  a  sound  speed 
of  c  =  1800  m/s.  The  canonical  Munk  SSP  from  Chapter  III.D  represents  the  sound 
speed  of  the  ocean  in  the  trench.  Sound  speed  increases  linearly  below  the  DSC 
corresponding  with  pressure  change  as  a  function  of  depth. 
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Figure  13.  Trench  A  bathymetry  representing  a  symmetrical  canyon. 


Bathymetry-  Subduction  Trench 


Figure  14.  Trench  B  bathymetry  representing  a  subduction  zone. 
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3. 


Results  for  Trench  A 


The  high  angled  walls  of  Trench  A  result  in  horizontal  surveillance  range 
perpendicular  to  the  trench  reduced  from  blockage  by  the  walls.  While  the  detection 
range  is  limited  by  this  effect,  the  walls  are  also  able  to  concentrate  detection  in  the  area 
that  is  still  visible.  Sound  that  does  not  reach  the  sensors  directly  at  lower  angles  can  still 
reach  the  sensor  by  reflection  from  the  angled  seafloor.  As  a  result,  the  sensor  may  see 
increased  sound  energy  from  the  target  at  certain  bearings  by  a  combination  of  direct  and 
bottom-reflected  paths.  Figure  15  shows  the  transmission  loss  for  a  flat  bottom  area 
while  Figure  16  shows  how  this  same  area  is  affected  by  the  trench.  Signal  drops  off 
outside  the  trench  due  to  blockage,  but  areas  near  the  walls  see  improved  signal  levels 
due  to  the  reflective  slope.  Figure  17  shows  that  for  a  target  at  10  m  below  the  surface, 
the  trench  provides  about  5  dB  in  improved  signal  energy  at  10  km  from  the  sensor. 
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Figure  15.  Constant  bathymetry  100  Hz  Transmission  Loss;  receiver  at  4390  m. 


o 

500 

1000 

1500 

2000 

2500 

3000 

3500 

4000 


BELLHOP-  Munk  Profile  Bottom  Source  -89  0.  f=100Hz,  Zs=4390m. 


01  15  20  25 

Range  (km) 


Figure  16.  Trench  A  100  Hz  Transmission  Loss;  receiver  at  4390  m. 
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BELLHOP-  Munk  Profile  Bottom  Source  -89  0:  f=100Hz.  Zs=4390m.  Zr=10m 


Figure  17.  Trench  A  100  Hz  Transmission  Loss;  source  at  10  m; 

receiver  at  4390  m. 


4.  Results  for  Trench  B 

The  lower  angle  walls  present  on  Trench  B  provide  less  focusing  of  the  sound 
energy  than  was  seen  with  Trench  A.  Figure  18  shows  that  the  gain  in  signal  at  distant 
ranges  is  less  significant.  With  the  steeper  angle  of  the  high  wall,  it  is  expected  that  its 
influence  would  be  more  significant  but,  as  Figure  19  shows,  there  is  still  little  effect. 

What  is  of  note,  in  Figure  19,  is  the  dramatic  drop  in  signal  level  toward  the 
trench  edge.  A  possible  explanation  is  that  the  shallower  depth  above  the  wall,  compared 
with  the  trench  bottom,  leaves  less  room  for  long-range  propagation  resulting  in  more 
scattering  and  signal  loss  before  reaching  the  sensor. 
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BELLHOP-  Munk  Profile  -89  0:  f=100Hz,  Zs=9140m,  Zr=10m 


Figure  18.  Trench  B  low  wall  100  Hz  Transmission  Loss;  source  at  10  m  depth; 

receiver  at  9140  m. 


BELLHOP-  Munk  Profile  -89  O:  f=100Hz,  Zs=9140m,  Zr=10m 


Figure  19.  Trench  A  high  wall  100  Hz  Transmission  Loss;  source  at  10  m  depth; 

receiver  at  9140  m. 
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5. 


Examples  and  Effects  of  Severe  Bathymetry 


Trenches  A  and  B  represent  likely  seafloor  geometries  to  be  found  in  many 
operating  environments.  In  some  locations,  it  may  be  possible  to  exploit  more  extreme 
geological  features  to  enhance  detection.  A  severe  trench  with  reflective  bathymetry  will 
limit  the  angle  of  coverage  to  the  shape  and  angle  of  the  trench  walls,  but  may  also 
greatly  focus  the  cone  of  detection  to  a  specific  area.  Trench  C,  represented  in  Figure  20, 
shows  an  example  of  a  narrow  trench  shaping  the  cone  of  detection.  Much  energy  is  lost 
to  reflection,  but  the  trench  does  provide  up  to  10  dB  of  improved  signal  strength  near  the 
vertical  axis. 


BELLHOP  100Hz  w  Bathymetry 


Figure  20.  Ray  trace  diagram  Trench  C  100  Hz;  receiver  at  3995  m. 
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BELLHOP-  Munk  Profile  Bottom  Source  -89  0.  f=100Hz,  Zs=3990m 
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Figure  21.  Trench  C  100  Hz  Transmission  Loss;  receiver  at  3995  m. 


BELLHOP-  Munk  Profile  Bottom  Source:  f=100Hz,  Zs=3990m,  Zr=10m 


Figure  22.  Trench  C  100  Hz  Transmission  Loss;  source  at  10  m;  receiver  at  3995  m. 
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More  significant  than  the  increase  in  signal,  is  the  reduction  in  ambient  noise 
expected  in  the  trench,  which  improves  SNR  by  decreasing  NL.  Figure  22  shows  a  large 
decrease  in  signal  energy  outside  the  trench  at  3  km  from  the  receiver.  This  narrower 
cone  of  detection  will  minimize  environmental  noise  received  and  will  not  see  contacts 
that  may  have  previously  been  visible  with  a  wider  cone.  This  can  help  in  isolating  and 
recognizing  a  target  of  interest  from  other  vessels  on  the  water. 

While  such  a  narrow  trench  can  provide  improved  SNR,  it  is  an  exceptional  case 
that  does  not  account  for  the  more  realistic  and  less  ideal  conditions  expected  in  an 
operational  environment.  A  trench  usually  will  not  be  symmetric  in  shape,  and  may 
likely  have  a  high  ridge  wall  opposite  a  relatively  open  ocean  environment,  as  seen  in 
Figure  23. 

For  modeling  Trench  D,  the  simulation  combines  the  cross-sections  of  a  low 
trench  (Figure  24)  with  the  steep  wall  of  the  Trench  C  (Figure  20).  The  high  wall  acts 
similar  to  the  narrow  trench  by  deflecting  noise  above  the  sensor,  leading  to  reduced  TL 
and  NL.  Figure  25  shows  improved  signal  strength  above  the  sensor  and  at  8500  m.  The 
lack  of  a  wall  opposite  the  high  ridge  prevents  reflection  towards  the  ridge  side.  The 
result  of  this  is  a  blind  spot  above  the  wall.  Figure  26  shows  that  detection  rapidly 
decreases  at  ranges  greater  than  2  km  from  the  sensor.  This  can  be  beneficial  in  contact 
separation  if  the  wall  is  blocking  a  shipping  lane  or  high-traffic  area  that  can  otherwise 
mask  target  detection. 
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Figure  23.  Ray  trace  diagram  Trench  D;  receiver  at  3995  m. 


BELLHOP-  Munk  Profile  Bottom  Source  Low  Trench  .  f=100Hz,  Zs=3990m. 
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Figure  24.  Trench  D  low  wall  100  Hz  Transmission  Loss;  receiver  at  3995  m. 
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Figure  25.  Trench  D  low  wall  100  Hz  Transmission  Loss;  source  at  10  m  depth; 

receiver  at  3995  m. 


BELLHOP-  Munk  Profile  Bottom  Source  Subduction  Trench:  f=100Hz,  Zs=3990m,  Zr=10m 


Figure  26.  Trench  D  high  wall  100  Hz  Transmission  Loss;  source  at  10  m  depth; 

receiver  at  3995  m. 
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V.  NOISE  EFFECT  ON  DETECTION 


A.  NOISE  DIRECTIONALITY  IN  THE  DEEP  OCEAN 

At  low  frequencies,  noise  arrives  at  a  bottom-mounted  deep  ocean  sensor 
predominantly  from  the  horizontal  direction  rather  than  from  overhead,  while  the 
opposite  is  true  at  high  frequencies.  Figure  27  shows  the  decibel  change  in  ambient  noise 
NL  as  a  result  of  vertical  beam  angle  (6)  at  the  receiver. 


(a)  (b) 


SPECTKUM  LEVEL  IN  UN,?  0»b)' /  STEH*DI*N  /  CYCLE  per  SECOND  SPECYBUM  LEVEL  IN  at//  (/lV]1  /  STERAOIAN  /  CYCLE  per  SECOND 

(*)  (0 

Figure  27.  Noise  Directionality  polar  plot  with  Beaufort  number  (from  [19]). 
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Noise  emanating  from  distant  sources  arrives  at  the  receiver  near  the  horizontal 
with  the  high  frequency  components  strongly  attenuated  in  transit.  Noise  from  the 
vertical  originates  from  sources  above  the  sensor  and  is  subject  to  much  less  attenuation 
by  virtue  of  the  shorter  distance.  Higher  frequency  noise  energy  arriving  at  near  vertical 
angles  apparently  originates  near  the  sea  surface  and  is  strongly  dependent  on  wind 
conditions,  as  measured  on  the  Beaufort  scale  [19]. 

B.  NOISE  CHARACTERISTICS  IN  A  TRENCH  ENVIRONMENT 

Since  propagation  of  lower  frequency  sound  is  sustained  at  great  distances  in  the 
ocean,  noise  energy  below  300  Hz  is  a  large  component  of  the  noise  received  at  the 
sensor.  Figure  6  shows  that  most  tonals  and  shipping  noise  are  at  low  frequencies. 
Therefore,  low-frequency  ambient  noise  can  adversely  affect  SNR.  Canyon  walls  and 
other  bathymetric  features  can  act  to  both  emphasize  a  signal  from  overhead  and,  at  the 
same  time,  suppress  noise  from  distant  sources. 

C.  CALCULATING  NOISE  LEVELS  AT  A  SENSOR 


Noise  spectrum  level  at  the  bottom-mounted  hydrophone  for  a  given  look  azimuth 
(f)  is  formulated  as: 


n_ 

N,=  ]  N(0)d0  (3) 

9M) 


where  6  is  the  vertical  look  angle,  0  -  0o  (</;)  is  the  elevation  slope  of  the  trench  wall  in  the 
look  azimuth  q J,  and  ()  =  Ky/1  is  the  vertical  upward  looking  direction. 

The  curve  in  Figure  28  is  converted  to  linear  intensity  from  decibels  (dB),  using 
Equation  4. 

NdB(e) 

NI(e)  =  10  10  (4) 


34 


Figure  28.  Graphical  representation  of  the  noise  calculation  showing  the  noise 

removed  by  the  trench  slope  0o. 


1.  Azimuth  Independent  Calculation 

The  simplest  assumption  for  calculating  the  noise  levels  in  a  trench  is  that  the 
sensor  is  blocked  on  all  sides  by  the  bathymetry.  This  is  similar  to  the  sensor  being 
placed  in  a  hole  in  the  ocean  floor  with  0o  constant  at  all  azimuths  (f).  Calculating  the 
noise  area  removed  by  this  method  simply  involves  treating  the  area  calculated  by 
Equation  3  as  a  solid  of  revolution  to  get  the  full  noise  volume. 

In  f 

Ntot  ~  f  f  N{d)  cos  0d8d(f)  (5) 

o  0M) 

where  Ntot  is  the  total  noise  volume  reaching  the  sensor  in  the  trench.  The  factor  cos  6 
accounts  for  the  spherical  nature  of  the  receiver  field  and  reduces  the  noise  concentration 
per  steradian  as  the  elevation  angle  increases. 
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2. 


Azimuth  Dependent  Calculation 


An  alternate  assumption  is  that  the  trench  is  a  parallel  trough  on  the  ocean  bottom 
and  the  area  of  noise  intensity  removed  depends  on  the  angle  (j)  to  the  sensor.  Therefore, 
N (0)  is  modeled  in  segments  of  n  intervals  over  a  quadrant  to  account  for  the  variations 
in  noise  volume. 


N,r 


In  2  2r 

—  ^  J  N(9)cos&d9 

2  n 


(6) 


n=number  of  segments 


Figure  29.  Overhead  graphical  representation  showing  the  segments  of  N,/,  at  each 
azimuth  (j)  which  make  up  the  total  volume  removed  by  the  trench  (Nu,t). 
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The  ratio  difference  between  Ntot  to  the  noise  level  without  the  trench  bathymetry 
(Nref)  represents  the  decibel  change  (AdB)  in  noise  caused  by  the  trench,  as  shown  in 
Equation  7. 


AdB  =  10  log 


\Nref  J 


(7) 


3.  Results 

The  azimuth-independent  noise  calculation  for  Trenches  A  and  B  shows  a 
decrease  in  noise  levels  of  4.6  dB  and  2.2  dB  respectively.  Since  low-frequency  noise 
arriving  from  horizontal  ranges  has  250  times  greater  intensity  than  overhead  sources, 
low  trench  walls  remove  a  significant  amount  of  the  total  noise  energy  from  the 
environment.  Trench  A  reduces  the  noise  level  by  65%  while  Trench  B  removes  40%  of 
the  noise  energy  from  the  environment. 

Using  the  azimuth-dependent  calculation  shows  a  decrease  of  2.4  dB  for  Trench 
A  and  1.4  dB  for  Trench  B.  Since  the  trench  slope  decreases  at  angles  near  (f>  =  90°  when 
using  the  azimuth  dependent  calculation,  it  is  expected  that  the  noise  reduction  will  not 
be  as  great  as  with  the  azimuth  independent  model. 

At  higher  frequencies  (>300  Hz),  where  the  maximum  noise  energy  begins  to 
arrive  from  steeper  vertical  angles,  it  is  not  expected  that  the  trench  will  cause  a 
significant  decrease  in  noise  energy. 

4.  Other  Examples 

Along  with  the  representative  trenches  above,  localized  canyons  and  ridges 
around  the  subduction  zones  have  the  potential  to  provide  steeper  angles  and  improved 
SNR.  Figure  30  shows  how  noise  level  changes  with  the  angle  from  the  sensor  to  the  top 
of  the  trench,  using  the  azimuth  dependent  calculation.  A  canyon  with  30-degree  walls 
can  decrease  noise  by  7.5  dB  and,  at  45  degrees,  the  reduction  is  10  dB.  Ridges,  where 
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there  is  only  a  single  wall  to  block  noise  (0iowerwaii  =  0°),  the  maximum  decrease  is  limited 
to  3  dB,  since  only  half  of  the  noise  field  can  be  blocked. 


Noise  Suppression  for  trench  placement  of  Sensor 


Figure  30.  Noise  suppression  data  for  sensor  placement  within  a  trench  for  symmetric 
(blue)  and  uneven  trenches  based  on  the  angle  of  the  lower  wall  (black). 


5.  Factors  Not  Considered  in  Noise  Calculations 

The  noise  mechanisms  in  Trench  B  (subduction  trench)  are  more  complex  than 

either  calculation  can  fully  account  for  due  to  its  vast  size  and  depth.  Pressure  is  the  main 

driver  of  sound  speed  below  the  DSC,  and  the  SSP  has  a  positive  speed  gradient  that 

deflects  sound  upward.  The  noise  calculation  for  both  models  assumes  isospeed  when 

calculating  directional  effects  of  noise  propagation.  However,  upward  deflecting  sound 

is  less  likely  to  be  blocked  by  low-slope  bathymetry,  meaning  that  the  noise  reduction 

calculation  for  Trench  B  is  likely  an  overestimation. 

38 


The  deep  placement  of  sensors  means  sound  emanating  from  near  the  surface 
needs  to  travel  nearly  vertical  to  reach  the  sensor  due  to  sound  speed  deflection.  As  such, 
noise  that  emanates  from  distant  sources  travels  on  a  path  that  undergoes  several  surface 
and  bottom  reflections  before  reaching  the  sensor.  Figure  19  shows  that  signal  level  in 
Trench  B  drops  significantly  on  the  continental  side  of  the  trench.  Shallow  areas  such  as 
coastal  water  will  subject  traveling  noise  to  greater  attenuation  than  deeper  ocean 
environments.  Therefore,  placing  the  sensor  low  in  a  trench  can  improve  ambient  noise 
levels  through  effects  other  than  blockage. 

Assuming  a  uniform  noise  environment,  the  high  wall  on  Trench  B  can  remove 
sound  from  half  the  field  giving  up  to  a  3dB  possible  reduction  in  noise.  However,  it  is 
likely  that  the  continental  shelf  will  have  higher  ambient  noise  levels  than  the  open  ocean 
due  to  a  greater  presence  of  local  traffic.  Therefore,  having  a  higher  wall  towards  the 
continental  shelf  means  the  stronger  contribution  to  the  ambient  noise  is  suppressed. 

Calculations  in  Matlab,  along  with  empirical  data  for  the  noise  curves  in  Figure 
27,  are  available  in  Appendices  A-D. 
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VI.  SENSOR-TO-SENSOR  COMMUNICATIONS 


A.  DEEP  SOUND  CHANNEL  PROPAGATION 

The  high  walls  and  ridge  formations  around  trenches  have  the  potential  to  create 
obstacles  to  modem-to-modem  communications  by  physical  obstruction  within  the  DSC. 
This  can  result  in  reflection,  attenuation,  and  scattering  of  the  signal,  reducing  the 
strength  and  hindering  communications. 

The  following  simulations  consider  a  5000  m  long  seamount  between  a  pair  of 
modems  obstructing  their  communication  channel.  The  seamount  is  10-15  km  from  the 
transmitting  modem.  The  maximum  depth  of  the  seamount  is  varied  to  assess  how  the 
overall  transmission  capability  is  affected. 

The  optimal  communications  modem  placement  is  at  or  just  below  the  DSC  axis. 
The  gateway  node  can  connect  either  to  this  DSC  node  or  directly  to  a  seabed  node.  The 
operating  frequencies  of  the  acoustic  modems  should  be  as  low  as  practical,  while  still 
retaining  adequate  spectral  bandwidth.  The  recommend  Deep  Seaweb  communications 
band  is  below  8  kHz.  Traditional  Seaweb  operating  bands  are  closer  to  12  kHz  [14]. 

1.  On-Axis  Modem  Results 

The  source  is  placed  on  the  DSC  axis  at  a  depth  of  1000  m,  corresponding  to  the 
optimum  location.  Figure  32  shows  the  different-sized  seamounts,  and  Figure  31  shows 
the  effect  each  has  on  the  communication  signal.  Obstructions  below  1000  m  do  not 
impact  the  signal  transmission  significantly.  This  is  mainly  due  to  the  fact  that,  at  this 
depth,  the  signal  can  travel  horizontally  through  the  DSC  channel  with  little  variation  in 
path.  The  interference  that  does  exist  comes  from  sporadic  destructive  interference  of  the 
signal  as  the  reflected  multipath  converges  with  the  non-reflected  path.  With  the 
seamount  at  shallower  depths  above  1000  m,  the  effects  become  more  apparent,  and  the 
obstructions  result  in  shadow  zones  where  no  signal  can  reach  the  receiver.  The  effect 
becomes  even  greater  when  the  seamount  cuts  through  a  larger  cross-section  of  the 
channel.  At  500  m,  the  channel  is  effectively  blocked. 
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BELLHOP-  Sound  Propogation  in  Channel:  f=10000Hz,  Zs=1000m,  Zr=1000m 
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Figure  3 1 .  Modem  10  kHz  TL  vs.  Range  for  obstructed  DSC  with  on-axis  modem. 


Bathymetry  Sound  Channel  Propogation 


Figure  32.  Bathymetry  vs.  Range  for  Figure  3 1 . 
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2. 


Offset- Axis  Modem  Results 


With  the  source  on  the  DSC  axis,  the  communication  channel  is  theoretically 
optimized  and  will  experience  the  least  transmission  loss.  This  placement  is  not  likely  in 
an  operational  environment  due  to  currents  that  will  cause  the  modem  to  shift  away  from 
the  DSC  axis.  It  is  known,  from  prior  testing,  that  placement  off  the  axis  results  in  a 
larger  refracted  path  with  increased  propagation  and  greater  attenuation  [14]. 
Propagation  is  modeled  for  a  source  at  1200  m  (offset  200  m  below  DSC  axis)  with 
seamount  heights  of  1500  m,  1000  m,  and  750  m.  Figure  33  shows  the  results  of  the  off- 
axis  simulation,  which  are  similar  to  the  on-axis  modeling  for  signal  blockage. 
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BELLHOP-  Sound  Propogation  in  Channel:  f=10000Hz,  Zs=1200m,  Zr=1000m 
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33.  Modem  10  kHz  TL  vs.  Range  for  obstructed  DSC  with  off-axis  modem. 


Bathymetry  Sound  Channel  for  Offset  Source 


Figure  34.  Bathymetry  vs.  Range  for  Figure  33. 
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VII.  CONCLUSION,  RECOMMENDATIONS,  WAY  AHEAD 


Seaweb  has  been  tested  extensively  in  various  stressful  environments.  The 
system  has  a  high  maturity  level  and  the  ability  to  detect  contacts  in  areas  with  high 
levels  of  background  noise.  Deep  Seaweb  is  a  more  recent  concept  that  uses  many  of  the 
advancements  of  Seaweb  technology  for  application  in  the  open  ocean  with  greater 
coverage.  Deep  Seaweb  exploits  greater  ocean  depths  and  sound  channel  features  for 
improved  sensor  performance  and  longer  communication  links. 

Exploiting  bathymetry  can  further  improve  Deep  Seaweb ’s  effectiveness. 
Simulation  results  show  that  placement  of  the  sensors  at  the  bottom  of  a  deep  ocean 
trench  improves  signal  reception  by  concentrating  and  amplifying  the  acoustic  signal 
from  a  vessel  of  interest.  Moreover,  trenches  provide  significantly  reduced  ambient  noise 
levels  by  blocking  horizontally  arriving  low-frequency  noise  resulting  in  improved  SNR 
for  better  detection  and  fewer  false  alarms. 

The  dynamics  and  environmental  characteristics  of  trenches  are  not  well 
understood.  Deep  placement  of  the  system  should  mitigate  variations  from  temperature 
and  pressure  on  the  propagation  of  sound  reaching  the  sensors.  The  noise  models  do  not 
account  for  flow  noise  or  seismic  events  possible  within  deep  ocean  trenches,  which  can 
increase  the  noise  level  at  the  sensor.  In  addition,  variations  in  the  bathymetry  or  the 
characteristics  of  the  ocean  floor  can  affect  the  performance  of  the  sensor  in  ways  not 
expected.  More  complex  trench  models  to  simulate  sound  and  noise  propagation  are 
necessary  to  better  understand  environmental  impact  on  the  system. 

Development  and  testing  of  Deep  Seaweb  in  an  operational  environment  is  the 
next  step  towards  deployment  of  the  system  as  an  autonomous  acoustic  network  that  can 
support  enforcement  of  maritime  law  and  interdiction  of  illegal  activities  in  open  waters. 
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APPENDIX  A.  AZIMUTH  INDEPENDENT  NOISE  MODEL 

MATLAB  CODE 


%This  program  calculates  the  change  in  noise  intensity  in 
%Trenches  A&B  using  azimuth  independence  and  a  -40dB 
%base  for  angles  blocked  by  the  trench  walls 

clear  all 

load  NoiseCurve 

%%  Define  trench  parameters 
L  1=1 3000;  %Trench  1  length  (m) 

Dl=2330;  %Trench  1  depth  (m) 

L2 1=48000;  %Trench  2  sidel  length  (m) 

L22=65000;  %Trench  2  side2  length  (m) 

D2 1=3350;  %Trench  2  side  1  depth  (m) 

D22=8850;  %Trench  2  side2  depth  (m) 

%%  Calculate  noise  from  graph 

%  file  'NoiseCurve'  contains  measurements  for  x  and  y  of  the  noise  curve.  Those 
%  values  are  referenced  to  values  of  dB 
db0=  1 .06;  %reference  length  for  dB=0 

db30=2. 12;  %reference  length  for  dB=30 

ddb=30/(db30-db0); 

x=(noisecurve( : ,  1 )); 
y=(noisecurve( :  ,2)); 

noisecurvel  12(:,1)  =  (0:.001:max(noisecurve(:,l))); 
noisecurvel  12(:, 2)=interpl(x,y, noisecurvel  12(:,1)); 

noisecurvel  12(:,3)=ddb.*noisecurvel  12(:,1); 
noisecurvel  12(:,4)=ddb.*noisecurvel  12(:,2); 

%  calculate  angle  phi  of  points  on  noise  curve 
phi=atan(noisecurvel  12(:,2)./noisecurvel  12(:,1)); 

%calculate  noise  intensity 

noisecurvel  12(:,5)=10.A(noisecurvel  12(:,3)-/10); 

noisecurvel  12(:,6)=10.A(noisecurvel  12(:,4)./10); 
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figure  (1);  plot  (noisecurvel  12(:,3),noisecurvel  12(:,4)); 
title  ('Reference  Noise  Level  (dB)') 

figure  (2);  plot  (noisecurvel  12(:, 5), noisecurvel  12(:,6)); 
title  ('Reference  Noise  Level  Intensity') 

%%  Defining  -40Db  power 

db40=0.000L*(max(noisecurvel  12(:,5)));  %noise  level  40dB  below  max  noise  level 
140=0. 5*pi*(db40)A2;  %40dB  area  of  a  half  disc 

%%  reference  noise 

Iref=pi.*max(noisecurvel  12(:,5)).*trapz(noisecurvel  12(:, 5), noisecurvel  12(:,6)); 

%%  Intergrating  Noise  Level  Trench  A 
n2=0; 
ltot=0; 

phi_t=atan(D  1/Ll); 
m=l; 

while  (phi(m,l)>phi_t) 
m=m+l; 
end 

phi_r=phi(m,l); 

I=pi.*noisecurvel  12(m,5).*(trapz(noisecurvel  12(l:m, 5), noisecurvel  12(l:m,6))- 
trapz(noisecurvel  12(1  :m,  5),  (noisecurvel  12(m,6)./noisecurvel  12(m,5)).*noisecurvel  12( 
l:m,5))); 

Itot=Itot+I; 

%Total  change  in  noise  level  in  dB 
dBchange=  1 0  *  log  1 0(Itot/Iref); 

disp( ['Noise  Level  decrease  for  Trench  A  ='  ,num2str(abs(dBchange)),  'dB']); 

%%  Integrating  Noise  Level  Trench  B 

%side  1 
n2=0; 

Itot2 1=0; 

phi_t=atan(D21/L2 1); 
m=l ; 

while  (phi(m,l)>phi_t) 
m=m+l; 
end 
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phi_r=phi(m,l); 

I=pi.*noisecurvel  12(m,5).*(trapz(noisecurvel  12(l:m,5),noisecurvel  12(l:m,6))- 
trapz(noisecurvel  12(1  :m,5),(noisecurvel  12(m,6)./noisecurvel  12(m,5)).*noisecurvel  12( 
1  :m,5)));+I40*phi_r/pi; 

Itot21=Itot21+I; 

%side  2 
n2=0; 

Itot22=0; 

phi_t=atan(D22/L22); 
m=l ; 

while  (phi(m,l)>phi_t) 
m=m+1 ; 
end 

phi_r=phi(m,l); 

I=pi.*noisecurvel  12(m,5).*(trapz(noisecurvel  12(l:m,5),noisecurvel  12(l:m,6))- 
trapz(noisecurvel  12(1  :m,5),(noisecurvel  12(m,6)./noisecurvel  12(m,5)).*noisecurvel  12( 
1  :m,5)))+I40*phi_r/pi; 

Itot22=Itot22+I; 

%Total  change  in  noise  level  in  dB 
dBchange=  1 0  *  log  1 0((Itot2 1  +Itot22)/ (2  *Iref)); 

disp( ['Noise  Level  decrease  for  Trench  B  ='  ,num2str(abs(dBchange)),  'dB']); 
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APPENDIX  B.  AZIMUTH  DEPENDENT  NOISE  MODEL  MATLAB 

CODE 


%This  program  calculates  the  change  in  noise  intensity  in 
%Trenches  A&B  using  azimuth  dependence  and  a  -40dB 
%base  for  angles  blocked  by  the  trench  walls 

clear  all 

load  NoiseCurve 


%%  Define  trench  parameters 
L  1=1 3000;  %Trench  1  length  (m) 

Dl=2330;  %Trench  1  depth  (m) 

L2 1=48000;  %Trench  2  sidel  length  (m) 

L22=65000;  %Trench  2  side2  length  (m) 

D2 1=3350;  %Trench  2  side  1  depth  (m) 

D22=8850;  %Trench  2  side2  depth  (m) 

%%  Calculate  noise  from  graph 

%  file  ’NoiseCurve’  contains  measurements  for  x  and  y  of  the  noise  curve.  Those 
%  values  are  referenced  to  values  of  dB/str 
db0=l  .06;  %reference  length  for  dB=0 
db30=2. 12;  %reference  length  for  dB=30 


ddb=30/(db30-db0); 


x=(noisecurve( ; ,  1 )); 
y=(noisecurve( ;  ,2)); 

noisecurvel  12(:,1)  =  (0:.001:max(noisecurve(;,l))); 
noisecurvel  12(:, 2)=interpl(x,y, noisecurvel  12(:,1)); 

noisecurvel  12(:,3)=ddb.*noisecurvel  12(;,1); 
noisecurvel  12(:,4)=ddb.*noisecurvel  12(:,2); 

%  calculate  angle  theta  of  points  on  noise  curve 
phi=atan(noisecurvel  12(:,2)./noisecurvel  12(;,1)); 

%calculate  noise  intensity 
noisecurvel  12(:,5)=10.A(noisecurvel  12(;,3)-/10); 
noisecurvel  12(:,6)=10.A(noisecurvel  12(:,4)./10); 
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figure  (1);  plot  (noisecurvel  12(:,3),noisecurvel  12(:,4)); 
title  ('Reference  Noise  Level  (dB)') 

figure  (2);  plot  (noisecurvel  12(:, 5), noisecurvel  12(:,6)); 
title  ('Reference  Noise  Level  Intensity') 


%%  Defining  -40dB  power 

db40=0.0001.*(max(noisecurvel  12(:,5)));  %noise  level  40dB  below  max  noise  level 
140=0.5  *pi*(db40)A2;  %40dB  area  of  a  half  disc 

%reference  noise 

Iref=pi.*max(noisecurvel  12(:, 5)). *trapz(noisecurvel  12(:, 5), noisecurvel  12(:,6)); 

%%  Integrating  Noise  Level  Trench  A 
n2=0; 
ltot=0; 

for  n=pi/180:pi/180:pi/2 
n2=n2+l; 

L=Ll/cos(n); 
phi_t=atan(D  1/L); 
m=l; 

while  (phi(m,l)>phi_t) 
m=m+l; 
end 

phi_r=phi(m,l); 

I=(pi.*noisecurvel  12(m,5).*(trapz(noisecurvel  12(1  :m, 5), noisecurvel  12(1  :m,  6))- 
trapz(noisecurvel  12(1  :m,  5),  (noisecurvel  12(m,6)./noisecurvel  12(m,5)).*noisecurvel  12 
(1  :m,5))))+I40*phi_r/pi; 

Itot=Itot+I; 

end 

%Calculating  total  noise  power 
Itot=Itot/n2; 

%Total  change  in  noise  level  in  dB 
dBchange=  1 0  *  log  1 0(Itot/Iref); 

disp([’Noise  Level  change  for  Trench  A  ='  ,num2str(abs(dBchange)),  'dB']); 

%%  Integrating  Noise  Level  Trench  B 
%side  1 
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n2=0; 

ltotl=0; 


for  n=pi/180:pi/180:pi/2 
n2=n2+l; 

L=L21/cos(n); 
phi_t=atan(D2 1/L); 
m=l ; 

while  (phi(m,l)>phi_t) 
m=m+l ; 
end 

phi_r=phi(m,l); 

I=(pi.*noisecurvel  12(m,5).*(trapz(noisecurvel  12(l:m,5),noisecurvel  12(1  :m, 6))- 
trapz(noisecurvel  12(1  :m,5),(noisecurvel  12(m,6)./noisecurvel  12(m,5)).*noisecurvel  12 
(1  :m,5))))+I40*phi_r/pi; 

Itotl=Itotl+I; 

end 

%Calculating  total  noise  power 
Itotl=Itotl/n2; 

%side  2 
n2=0; 

Itot2=0; 

for  n=pi/180:pi/180:pi/2 
n2=n2+l; 

L=L22/cos(n); 
phi_t=atan(D22/L); 
m=l ; 

while  (phi(m,l)>phi_t) 
m=m+l ; 
end 

phi_r=phi(m,l); 

I=(pi.*noisecurvel  12(m,5).*(trapz(noisecurvel  12(l:m,5),noisecurvel  12(l:m,6))- 
trapz(noisecurvel  12(1  :m,5),(noisecurvel  12(m,6)./noisecurvel  12(m,5)).*noisecurvel  12 
(1  :m,5))))+I40*phi_r/pi; 


55 


Itot2=Itot2+I; 

end 

%Calculating  total  noise  power 
Itot2=Itot2/n2; 

Itot=(Itot  1  +Itot2)/2 ; 

%Total  change  in  noise  level  in  dB 
dBchange=  1 0  *  log  1 0(Itot/Iref); 

disp(['Noise  Level  change  for  Trench  B  ='  ,num2str(abs(dBchange)),  ’dB’]); 
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APPENDIX  C.  NOISE  SUPPRESION  GRAPH  MATLAB  CODE 


%This  program  calculates  the  change  in  noise  energy  level  in 
%trenches  over  a  range  of  angles  and  creates  a  graph  of  the  numbers 
%using  a  -40dB  base  for  angles  blocked  by  the  trench  walls 

clear  all 

load  NoiseCurve 
%%  Calculate  noise  from  graph 

%  file  ’NoiseCurve'  contains  measurements  for  x  and  y  of  the  noise  curve.  Those 
%  values  are  referenced  to  values  of  dB 
dbO=l  .06;  %reference  length  for  dB=0 
db30=2. 12;  %reference  length  for  dB=30 

ddb=30/(db30-db0); 

x=(noisecurve( : ,  1 )); 
y=(noisecurve( :  ,2)); 

noisecurvel  12(:,1)  =  (0:.001:max(noisecurve(:,l))); 
noisecurvel  12(:, 2)=interpl(x,y, noisecurvel  12(:,1)); 

noisecurvel  12(:3)=ddb.*noisecurvel  12(:,1); 
noisecurvel  12(:,4)=ddb.*noisecurvel  12(:,2); 

%  calculate  angle  phi  of  points  on  noise  curve 
phi=atan(noisecurvel  12(:,2)./noisecurvel  12(:,1)); 

%calculate  noise  intensity 
noisecurvel  12(:,5)=10.A(noisecurvel  12(:,3)-/10); 
noisecurvel  12(:,6)=10.A(noisecurvel  12(:,4)./10); 

%%  Defining  -40Db  power 

db40=0.000L*(max(noisecurvel  12(:,5)));  %noise  level  40dB  below  max  noise  level 
140=0. 5*pi*(db40)A2;  %40dB  area  of  a  half  disc 

%%  Integrating  Noise  Level  Trench  A 

/Preference  noise 

Iref=pi.*max(noisecurvel  12(:,  5)).  *trapz(noisecurvel  12(:,  5),  noisecurvel  12(:,6)); 
k=0; 
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for  phi_i=0:pi/36: 10*pi/36 
k=k+l; 
n=0; 
ltot=0; 

for  theta=pi/1800:pi/1800:pi/2 
n=n+ 1 ; 

phi_t=atan(tan(phi_i)*cos(theta)); 
m=l ; 

while  (phi(m,l)>phi_t) 
m=m+ 1 ; 
end 

phi_r=phi(m,l); 

I=(pi.*noisecurvel  12(m,5).*(trapz(noisecurvel  12(l:m,5),noisecurvel  12(l:m,6))- 
trapz(noisecurvel  12(1  :m,5),(noisecurvel  12(m,6)./noisecurvel  12(m,5)).*noisecurvel  12 
(1  :m,5))))+I40*phi_r/pi; 

Itot=Itot+I; 

end 

%Calculating  total  noise  intensity 
Itot=Itot/n; 

Itot_t(k,l)=Itot; 

Itot_t(l,l)=Iref; 

%%  Build  Table  for  Values  of  DB  vs.  Phi 

dBchange(k,  l)=phi_i*  1 80/pi;  %phi  values 
dBchange(k,2)=abs(10*logl0(Itot/Iref));  %dB  values 

j=i; 

while  j<=k 

Itot_t(j  ,k)=(Itot_t(k,  1  )+Itot_t(j ,  1  ))/2 ; 
dBchange(j,k+2)=abs(  1 0*log  1 0(ltot_t(j,k)/Iref)); 

j=j+i; 

end 

end 

%%  Build  Plot 

figure  (3);  plot(dBchange(: ,2),dBchange(: ,  1  ),'-bo ','LineWidth',3 •  • 
’MarkerEdgeColorVk','MarkerSize',5); 

title('Noise  Suppression  for  Trench  Placement  of  Sensor','fontsize',12); 
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xlabel('abs(NL_f-NL_o)  (dB) ’,'fontsize',  1 0); 

ylabel(' Approach  angle  of  bathymetry  wrt  sensor  (\phi)',Tontsize',10); 


hold  all 
for  d=l:l:k 

plot  (dBchange(l :d,d+2),dBchange(l :d,l),’-kVlinewidth',2); 
end 

grid  on; 
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APPENDIX  D.  NOISECURVE  DATA 


Calculations  for  noise  levels  using  code  in  Appendix  A-C  is  based  on  the  1 12Hz 
graph  in  Figure  27  for  Beaufort  <4  wind  levels.  The  following  table  is  a  representation  of 
this  graph  in  x-y  coordinates  and  loaded  in  file  Noisecurve.m  utilized  by  those  codes. 


X 

9 

x  (db) 

y  (db) 

1.95 

0 

55.18867925 

0 

1.94 

0.05 

54.90566038 

1.41509434 

1.91 

0.1 

54.05660377 

2.830188679 

1.9 

0.15 

53.77358491 

4.245283019 

1.88 

0.2 

53.20754717 

5.660377358 

1.86 

0.25 

52.64150943 

7.075471698 

1.83 

0.3 

51.79245283 

8.490566038 

1.8 

0.35 

50.94339623 

9.905660377 

1.76 

0.41 

49.81132075 

11.60377358 

1.72 

0.46 

48.67924528 

13.01886792 

1.68 

0.51 

47.54716981 

14.43396226 

1.63 

0.57 

46.13207547 

16.13207547 

1.58 

0.63 

44.71698113 

17.83018868 

1.54 

0.67 

43.58490566 

18.96226415 

1.48 

0.71 

41.88679245 

20.09433962 

1.43 

0.75 

40.47169811 

21.22641509 

1.39 

0.79 

39.33962264 

22.35849057 

1.32 

0.83 

37.35849057 

23.49056604 

1.26 

0.86 

35.66037736 

24.33962264 

1.21 

0.89 

34.24528302 

25.18867925 

1.15 

0.92 

32.54716981 

26.03773585 

1.08 

0.94 

30.56603774 

26.60377358 

1.01 

0.97 

28.58490566 

27.45283019 

0.94 

1 

26.60377358 

28.30188679 

0.86 

1.02 

24.33962264 

28.86792453 

0.77 

1.05 

21.79245283 

29.71698113 

0.68 

1.07 

19.24528302 

30.28301887 

0.61 

1.08 

17.26415094 

30.56603774 

0.51 

1.09 

14.43396226 

30.8490566 

0.45 

1.095 

12.73584906 

30.99056604 

0.38 

1.095 

10.75471698 

30.99056604 

0.33 

1.1 

9.339622642 

31.13207547 

0.28 

1.1 

7.924528302 

31.13207547 

0.23 

1.1 

6.509433962 

31.13207547 

0.18 

1.1 

5.094339623 

31.13207547 

0.13 

1.1 

3.679245283 

31.13207547 

0.08 

1.1 

2.264150943 

31.13207547 

0.03 

1.1 

0.849056604 

31.13207547 

0 

1.1 

0 

31.13207547 
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